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Abstract: We use the recently reported KamLAND measurements on oscillations 
of reactor anti-neutrinos, together with the data of previously reported solar neu- 
trino experiments, to show that: (1) the total 8 B neutrino flux emitted by the Sun is 
1.00(1 ±0.06) (lor) of the standard solar model (BP00) predicted flux, (2) the Kam- 
LAND measurements reduce the area of the globally allowed oscillation regions that 
must be explored in model fitting by six orders of magnitude in the Am 2 — tan 2 9 
plane, (3) LMA is now the unique oscillation solution to a CL of 4.7cr, (4) max- 
imal mixing is disfavored at 3.1a, (5) active-sterile admixtures are constrained to 
sin 2 7/ < 0.13 at la, (6) the observed 8 B flux that is in the form of sterile neutri- 
nos is O.OOIq^q (lo")) of the standard solar model (BP00) predicted flux, and (7) 
non-standard solar models that were invented to avoid completely solar neutrino 
oscillations are excluded by KamLAND plus solar data at 7.9a. We also refine quan- 
titative predictions for future 7 Be and p — p solar neutrino experiments. 
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1. Introduction 

KamLAND is the first experiment to explore with a terrestrial beam the region of 
neutrino oscillation parameters that is relevant for solar neutrino oscillations. This 
spectacular achievement [|TJ makes possible new studies of the physics of neutrinos 
[0, [| [|, [| |7], [| [TU[ and of the solar interior [O] . We concentrate here on refining 
the allowed regions in neutrino oscillation space and on what can be learned about 
the total 8 B solar neutrino flux, as well as the sterile neutrino component of the 8 B 
neutrino flux. The procedures we use have been described in refs. [§, [IT2] . 



We assume throughout this paper that CPT is satisfied and that the anti- 
neutrino measurements by KamLAND apply directly to the neutrino section. In 
fact, the first KamLAND results already show that CPT is satisfied in the weak 
sector to a characteristic accuracy of about 10 -20 GeV [ffl. 
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We begin in section |2] by determining the allowed oscillation regions using the re- 
cently released KamLAND data. Our principal results, which are in good agreement 
with the analysis of the KamLAND collaboration |IJ, are summarized in figure [j]and 
in table |3|. We also show that the numerical values of the physical quantities that we 
calculate in this paper are insensitive to the details of the analysis of the KamLAND 
data (see discussion at the end of section |6j). 

We derive in section || a global solution for the neutrino oscillation parameters 
using, together with the KamLAND measurements, all the available solar neutrino 



data from the chlorine [13], gallium [14, 15,[L6[, Super-Kamiokande [17|, and SNO [18 



191 experiments. The enormous reduction in the size of the allowed oscillation regions 



is evident from a comparison of the 'Before' and 'After' allowed regions, figure |]and 
figure |j. 

We determine in section [| the total 8 B solar neutrino flux, as well as the uncer- 
tainties in the flux that are implied by the existing experimental data. This flux can 
be compared directly with the flux predicted by the standard solar model. 

Could there be a large but unobservable flux of 8 B solar neutrinos that reaches 
Earth in the form of sterile (right-handed) neutrinos? Prior to the announcement of 
the KamLAND results U, the answer to this question was a resounding "Yes" @, |(] 



At 3<T, the sterile component of the 8 B flux could be, prior to KamLAND, as large 
as 1.1 times the BPOO predicted flux, which amounts to 50% of the total flux for this 
extreme case. However, the KamLAND measurement allows us to place a strong 
upper limit on the flux of sterile 8 B neutrinos. We determine in section |5| the best- 
estimate sterile 8 B neutrino flux and the associated uncertainties in the sterile flux. 
We summarize and discuss our main conclusions in section 161. 



2. KamLAND only 

Figure [l] shows the allowed regions for anti-neutrino oscillations that we have de- 
termined using the description of the detector operation and the measured data 
presented in the recent KamLAND publication [|TJ. The left hand panel of figure [I] 
shows the allowed regions calculated using only the KamLAND total event rate; 
the right hand panel shows the allowed regions calculated using the experimental 
energy spectrum, which also contains a rate normalization. In computing the right 
hand panel of figure [TJ, we also used the data from the CHOOZ reactor experiment. 
This experiment rules out solutions with Am 2 greater than ~ 1 x 10~ 3 eV 2 at 3cr or 
~ 8 x l(T 4 eV 2 at 99% CL. 

Section [2.1| contains our version (used in the following sections of the paper) 
of the science results from the KamLAND measurements and may be of general 
interest. Section |2.2j contains a description of how we do the data analysis and is 
rather specialized. We describe in sec. |2.3| Earth matter effects in the KamLAND 



experiment. The results of sec. [2.3| show that matter effects must be included in 
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Figure 1: Anti-neutrino oscillation parameters regions allowed by KamLAND. 

The figure shows the allowed regions for anti-neutrino oscillation, v e -/-> P e , as implied by 
the first KamLAND results . The contours shown in the left hand panel were calculated 
using only the total rate observed in the KamLAND experiment; the contours in the right 
hand panel were calculated using the measured KamLAND energy spectrum. The contours 
shown in figure [l] are 90%, 95%, 99%, and 99.73% (3<r)confidence limits. The dark crosses 
mark the best-fit values for the energy spectrum. The non-observation of anti-neutrino 
oscillations in the CHOOZ |21] reactor experiment rules out solutions with Am 2 greater 
than ~ 1 x l(T 3 eV 2 at 3a or ~ 8 x l(T 4 eV 2 at 99% CL . 



future precision analyses of the KamLAND experiment. We have carried out the 
KamLAND-only and the solar plus KamLAND analyses described in this paper 
both with and without including matter effects for the KamLAND experiment and 
have verified that, at the present level of accuracy for KamLAND measurements, the 
predictions for future experimental measurements are unaffected by including matter 



effects for KamLAND. Section [2.2| and Section [2.3| are rather technical and can be 
skipped by the reader who is primarily interested in the science. 

We have verified that different plausible approaches to analyzing the KamLAND 
data do not lead to significant differences in the quantities in which we are most 
interested: the total 8 B solar neutrino flux, the sterile component of the 8 B neutrino 
flux, the CL at which maximal mixing is excluded, and the predictions for future 
solar neutrino experiments. 

2.1 The KamLAND allowed regions 

As expected, the allowed regions shown in figure ffl are in agreement with those 
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reported by the KamLAND collaboration. It was necessary to establish explicitly 
this agreement in order to have confidence in the subsequent steps of our analysis, 
which involve determining: the globally allowed regions (section |3.1.1|) , predictions 
for future solar neutrino experiments (section |3.1.2| ), the total 8 B solar neutrino flux 



(section ^|), an upper limit to the sterile component of the 8 B neutrino flux (section |5|), 
and an upper limit to the CNO contribution of the solar luminosity (ref. ||11|| ). 

The results we obtain for the global solar plus KamLAND contours of the al- 
lowed regions, and all of the quantities derived from these contours, are essentially 
independent of the details of the analysis of the KamLAND results. We have ver- 
ified this independence by analyzing the KamLAND measurements in a number of 
different ways (see discussion at the end of section 6). 

The most important aspect of figure [I] is the demonstration by the KamLAND 
results that anti-neutrinos oscillate with parameters that are consistent with the 
LMA solar neutrino solution ||22|| . We shall exploit this extraordinary achievement 
in detail in the following sections. The symmetry shown in figure [I] with respect to 
the line 9 = 7r/4 (or tan 2 9 = 1) is just a reflection of the fact that the KamLAND 
experiment is an essentially vacuum experiment (but cf. sec. |2]^ and fig. |^) and 



hence the survival probabilities depend almost entirely upon sin 2 29. 

For the rates-only analysis, the best-fit oscillation parameters form a degenerate 
line in the Am 2 -tan 2 # plane. For large Am 2 , the line only depends upon 9 and is 
described by tan 2 9 = 0.36. 

Expressed in terms of the rate expected if there were no flavor changes, the 
KamLAND rate is 

^KamLAND = 0.611 ± 0.094 . (2.1) 

The measured rate is in excellent agreement with the pre-KamLAND expectations 
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-^KamLAND, expected — 0-58_ ' 27 , (2.2) 

based on purely solar neutrino experiments. 

As can be seen in the right hand panel of figure [I] and as is summarized in 
table [1], the KamLAND spectral data lead to three local minima. The allowed region 
is separated into 'islands.' These islands correspond to oscillations with wavelengths 
that are approximately tuned to the average distance between the reactors and the 
detector, 180 km. The best-fit point in the lowest mass island (Am 2 = 1.5 x 10~ 5 eV 2 ) 
is near the first-maximum in the oscillation probability (minimum in the event rate). 
The overall best-fit point (Am 2 = 7.1 x 10 _5 eV 2 ) lies within the island around the 
second maximum in the oscillation probability. The tuning is less accurate for the 
best-fit point in the highest mass island. 

For each of the allowed islands in neutrino oscillation parameter space that are 
shown in the right hand panel of figure [I] (analysis of the energy spectrum), table [I] 
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Am 2 


tan 2 6 


Ay 2 • 

Arm n 


7.1 x 1(T 5 


5.2 x HT 1 (1.9 x 10°) 


0.0 


1.7 x 1(T 4 


3.5 x KT 1 (2.8 x 10°) 


1.6 


1.5 x 1(T 5 


3.7 x HT 1 (2.7 x 10°) 


3.0 



Table 1: KamLAND allowed regions. The table presents the best-fit oscillation pa- 
rameters for the current 3<r-allowed LMA islands determined by the KamLAND energy 
spectrum (see figure [l]). We also give the x 2 difference between the local minimum in each 
region and the global minimum. 



gives the best-fit values of Am 2 and tan 2 9. The global best-fit point for the Kam- 
LAND data set lies at Am 2 = 7.1 x 10~ 5 eV 2 and mixing angle tan 2 6 = 0.52. For 
the same mass difference, maximal mixing is slightly less favored, A% 2 = 0.4. With 
the present statistical accuracy, KamLAND cannot discriminate between large angle 
(non-maximal) mixing and maximal mixing. 

At what level are other, non-LMA oscillation solutions excluded? For Am 2 < 
10 _6 eV 2 , there is no significant oscillation probability in the KamLAND experiment. 
All such low Am 2 solutions (like the previously-discussed LOW and Vacuum solar 
neutrino oscillation solutions), have A% 2 > 15.9 and are excluded by at least 3.6a (2 
dof). 

Comparing the total number of observed and expected events in KamLAND 
we find that alternative explanations to the solar neutrino problem that predict no 
deficit in the KamLAND experiment are now excluded at 3.6a\ This applies, for 
instance, to spin flavor precession, or flavor changing neutrino interactions which are 
now ruled out at 3.6<r as the dominant mechanism for solar neutrino flavor conver- 
sion. In particular, all 'non-standard solar model' explanations of the 'solar neutrino 
problem' which assume particle physics with P ee = 1 (no neutrino oscillations) are 
also excluded by the KamLAND experiment at 3.6<r. 

2.2 KamLAND analysis 

^From the experimental point of view, KamLAND involves different systematic un- 
certainties than solar neutrino experiments. From a theoretical point of view, Kam- 
LAND uses anti-neutrinos to perform a pure disappearance experiment in vacuum. 
Solar neutrino studies use neutrinos to perform disappearance (or appearance) ex- 
periments in matter (or in vacuum). The analysis of the KamLAND data requires 
a treatment that is different from, and in many ways more simple than, the global 
analysis of solar neutrino data. 

Using the data provided in the KamLAND paper M, we calculate the positron 



spectrum in KamLAND detector with the procedures described in refs. B 53, E4L [25 



26] . In the absence of neutrino oscillations, we find (in agreement with ref. [|1|]) 86.8 
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expected neutrino events above 2.6 MeV visible energy, for 145.1 days of live time, 
a 408 ton fiducial mass with 3.46 xlO 31 free target protons, detection efficiency of 
78.3%, energy resolution a(E)/E = 7.5%/ 'yE, integrated total thermal power flux 
during the measurement live time of 254 Joule/cm 2 , and relative fission yields from 
different fuel components as specified in ref. ]TJ. The positron energy spectrum that 
we calculate is in excellent agreement with the energy spectrum without oscillations 
presented by the KamLAND collaboration. 

We perform, following the KamLAND collaboration, two analyses of the first 
KamLAND results: 1) including only the total event rate, and 2) including the total 
measured energy spectrum. For the rate only analysis, the \ 2 can be written as JD], 



2 (0.611 -R(6, Am 2 )) 2 

XRate 



a 



(2.3) 



R 



with <t^ 



(0.085) 2 + (0.041) 2 and 



, „ , n / Am 2 !/ 
R=l — sin 2 26 sin 2 



) 



;2.4) 



We denote by () the convolution over the energy spectrum, the interaction cross 
section, and the energy resolution function. For simplicity in the notation, we have 
neglected in writing eq. |2]^ all matter effects (see sec. |2.3| for details of the matter 
effects) . 

In the left hand panel of figure |IJ, we show the allowed regions (2 dof ) that are 
obtained from this rate-only analysis. The vertical asymptotes can be calculated 
simply from eq. (|2.3| ) since for large Am 2 the quantity ^sin 2 (Am 2 L/AE u )^ averages 
to 0.5 . Thus in the large Am 2 limit, the mixing angle boundary of the allowed region 
satisfies sin 2 (26 l ) = — 2cr R y / Ax 2 7L + 0.778. Here A% 2 is fixed by the CL of interest; 
Ax 2 = 5.99(11.83) for 95% CL (for 3a CL). 

We included the data from the KamLAND energy spectrum in a way that has 
become conventional when analyzing low statistics binned data |[27|| , assuming that 
data is Poisson distributed. We define the \ 2 function as 



X 2 (Am 2 



mm, 



E 



-2£ 



2(aiVf (Am 2 , i 



n: 



exp\ 



+ 2n; 



exp 




[a 



a, 



(2.5) 



MVS 



where a is an absolute normalization constant and cr sys = 6.42% represents the 
systematic uncertainties presented in table 2 of ref. 0. The KamLAND collaboration 
analyzed the event sample by separating the total rate from the spectrum shape and 
using a maximum likelihood method. 

In the right hand panel of figure [1], we show the allowed regions (2 dof) that we 
obtain from our analysis of the KamLAND spectrum. 



6 



We have performed a series of tests of the dependence of the allowed solutions on 
the statistical treatment of the data. For example, we have performed a standard x 2 
analysis using the first 5 energy bins above 2.6 MeV visible energy as presented by 
the KamLAND collaboration, but with all the events above the fifth bin (E > 4.725 
MeV) collected into a single bin (or into two bins). We then assumed Gaussian 
statistics for all bins. We find that the exact size and shape of the 3a allowed regions 
of the KamLAND-only analysis depend mildly on the strategy adopted, i.e., whether 
we use eq. (|2.5|) or the standard x 2 . However, the differences are unimportant for 
quantities calculated using the global solar plus KamLAND allowed regions. We 
have also studied the effect of carefully treating the energy dependence of the energy 
threshold error by computing, as a function of the oscillation parameters, the sys- 
tematic uncertainty for each bin due to the estimated systematic error for the energy 
scale (1.91% at 2.6 MeV(l|]). This refinement has no discernible effect on any of the 
quantities we calculate in this paper. 

We discuss at the end of section |6| the small quantitative dependence upon anal- 
ysis strategy of each of the quantities we calculate in this paper. 

2.3 Matter effects in the KamLAND experiment 

It is conventional to ignore matter effects in analyzing the KamLAND results. How- 
ever, as we shall see below, matter effects can influence the calculated anti-neutrino 
event rate by as much as 4% if the mass and mixing angle are near the region of 
maximum sensitivity. Thus matter effects will have to be included in future analy- 
ses of the KamLAND results when greater statistical and systematic precision are 
obtained. 

In the results described in this paper, we have included matter effects by assum- 
ing that the anti-neutrinos traverse a constant density medium (p = 2.7gmcm~ 3 , 
cf. ref. [^])with the path length appropriate to each reactor. Since matter effects 
are small, the assumption of a constant density is an excellent approximation. Once 
matter effects are included, the survival probability for the KamLAND experiment 
also depends on the active-sterile admixture. However, allowing sterile neutrinos 
reduces the influence of matter effects with respect to what is calculated for the pure 
active case. 

To the accuracy that can be discerned by even the most expert eye, matter effects 
do not change the contours of figures given in the present paper. 

We define, F(matter vs vacuum), as the fractional difference in the event rate 
for the KamLAND detector calculated with and without including matter effects in 
the Earth. Thus 

Ft \ -^matter -^vacuum «\ 
(matter vs vacuum) = . (2.6) 

R vacuum 

Figure |2] shows the contours of F(matter vs vacuum) in the Am 2 -tan 2 # plane 
for active neutrinos. The maximum value of |F(matter vs vacuum) | corresponds to 
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Figure 2: Matter effects in the KamLAND reactor experiment. The contours 
depict the fractional difference, defined in eq. 2.6, in the event rate calculated with and 
without including matter effects from anti-neutrinos traversing the Earth. 



-4.2% (+4.3%) at Am 2 = 3.55 x KT 5 eV 2 and tan 2 # = 0.43(2.18). Within the 
la (3a) allowed region of the analysis of the KamLAND energy spectrum (cf. the 
right hand panel of fig. [I]), the maximum value of |F(matter vs vacuum) | corresponds 
to 0.9% (3.9%). The maximum change in x 2 due to including matter effects in the 
KamLAND-only analysis is 0.2 (1.5) at la (3a). At the best-fit point for the analysis 
of the KamLAND energy spectrum, matter effects change the value of x 2 by only 
0.1%. 

We anticipate the results given latter in sec. ^| (see especially fig. f|) by noting that 
the matter effects in the KamLAND experiment are somewhat less important within 
the globally allowed solar plus KamLAND oscillation regions. For the global solution, 
|F(matter vs vacuum) | varies between 0.9% (2.7%) at la(3a). The maximum change 
in x 2 m the global analysis due to matter effects is 0.2 (1.0). 

For simplicity, we have described in the rest of this paper the KamLAND exper- 
iment as if matter effects could be completely ignored. 
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3. Global solutions: solar plus KamLAND 



We present and discuss in section |3J] the allowed regions in neutrino oscillation 
parameter space that are permitted at 3<r by the combined solar and KamLAND 
experimental data (see especially figure |]). We also present in this section (see espe- 
cially table the predictions of future solar neutrino observables that are calculated 
using the current global oscillation solutions. 

In order to isolate the crucial measurements that have identified LMA as the 



correct description of solar neutrino oscillations, we present in section |3.2| the results 
of an analysis of just the rates of the solar neutrino experiments. We show in this 
section that, if we did not have the essentially undistorted electron recoil energy spec- 
trum measured by Super-Kamiokande and confirmed by SNO, vacuum oscillations 
would be the preferred solution prior to the KamLAND measurements. 

We include two sections that are primarily intended for aficionados of neutrino 
oscillation analyses. We outline briefly in section |3.3| how we have combined the 



solar and the KamLAND analyses. We have already described in section 2.2 how 



we have analyzed the KamLAND measurements. In section |3.3.1 



we summarize our 



treatment of the solar data. We describe in section |3.3.2| our treatment of the sterile 
component of the solar neutrino flux. 

3.1 Scientific results 



We describe in section |3.1.1| the allowed regions in neutrino oscillation space that are 



determined by the solar and KamLAND data and in section |3.1.2j we present predic 



tions for future solar neutrino experiments of the presently allowed global oscillation 
solutions. 

3.1.1 Allowed regions 

Figure ^ presents the 'Before' figure, the allowed neutrino oscillation regions for all 
the available solar neutrino data prior to the announcement of the first KamLAND 
results. We show in figure |3] the results of our analysis for this paper of the pre- 
KamLAND solar neutrino data. We include here for the first time the recently 
reported GNO measurements [M; this inclusion does not affect figure 131 significantly. 



Our calculational procedures, which are described in ref. |12||, contain some refine- 



ments not taken account of in other analyses. However, the overall results for the 



pre-KamLAND allowed regions are very similar in most treatments, cf. [[28[ . 

Figure £| displays the 'After' figure, the much reduced allowed regions that are de- 
termined by including the recently announced KamLAND measurements [0 together 
with the solar neutrino data. 

Comparing figure [| with figure |], we see that KamLAND has shrunk enormously 
the allowed parameter space that is accessible to theoretical models, by five orders 
of magnitude in Am 2 and one order of mae nitude in tan 2 9. The LMA region m 
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Figure 3: Before: Solar neutrino oscillations before KamLAND. The contours 
shown in figure || are 90%, 95%, 99%, and 99.73% (3a")confidence limits. The global best- 
fit points are marked by a star. The input data used in constructing the figure include all 
the available solar neutrino measurements [1^, 18, 19, 14, 16|, 17] and the neutrino fluxes 
and uncertainties predicted by the BPOO solar model (2t| for all but the 8 B flux. The 8 B 
neutrino flux is treated as a free parameter. 



parameter space is the only remaining allowed solution. The best-fit LOW solution 
is excluded at 4.8a and vacuum solutions are excluded at 4.9a. The once 'favored' 
SMA solution is now excluded at 6.1a. 

The allowed range for the mixing angle within the la (3a) region shown in the 
'After' figure f| is 

0.36( 0.27) < tan 2 9 < 0.58 (0.94) . (3.1) 

For the pre-KamLAND analysis, figure |^, the allowed 3a domain for the mixing angle 
was 0.26 < tan 2 9 < 0.98. Therefore, the allowed range of mixing angles has only 
been marginally improved by the inclusion of the KamLAND data. 

At what confidence level is maximal mixing excluded? There are no solutions 
with tan 2 6» = 1 at the 3.1cr CL. For the pre-KamLAND analysis, figure |3|, there 
were no solutions with maximal mixing within the LMA region at 3.1a. However, 
maximal mixing was still allowed at 2.8a within the vacuum solution regions, which 
is now excluded. 

Figure |^ shows how solar neutrino experiments break the symmetry of the sur- 
vival probability, P cc (E), about a mixing angle of 9 = tc/2 that exists for vacuum 
oscillation experiments. Vacuum oscillation survival probabilities are identical for 
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Figure 4: After: Solar neutrino oscillations after KamLAND. The allowed con- 
tours shown in figure || include the recent KamLAND measurements and all the available 
solar neutrino data [|^, 14, 15, |1| [l?], [l]|, 19]. The contour limits are the same as in figure [j] 
and figure [|. 

two mixing angles 9i and 8 2 , with 6 2 = 7r/2 — 9i and a fixed value for Am 2 . The left 
hand panel shows the vacuum survival probabilities for a source-detector separation 
of 180 km, similar to the actual experimental situation for KamLAND. The right 
hand panels displays the daytime survival probabilities for the two LMA solar neu- 
trino oscillation solutions. Matter effects in the Sun are responsible for the different 
LMA survival probabilities. 

Which solar neutrino experiments choose the lower probabilities (the smaller 
value for tan 2 #)? Originally, the chlorine experiment, in conjunction with the stan- 
dard solar model predictions, indicated strongly that the smaller probabilities were 
correct. The measured event rate in the chlorine experiment, ~ 2.6 SNU was 
closer to 0.35 of the standard solar model prediction than to 0.65 of the solar model 
prediction (~ 8 SNU) p9| . This result was confirmed by the comparison of the Super- 
Kamiokande v — e scattering measurements with the SNO CC measurement [|3(J. The 
comparison within the same experiment of the SNO CC and NC rates selects unam- 
biguously the curve with the lower survival probabilities [18| . 

Finally, we have evaluated the present status of 'non-standard' solar model ex- 
planations of the solar plus KamLAND data under the hypothesis of standard par- 



ticle physics. We perform a fit to the total measured rates in the chlorine [ITS 
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E (MeV) 



Figure 5: Breaking the degeneracy: P cc vs E. Figure [| shows the survival probability, 
P e ej as a function of neutrino energy, E, for two mixing angles related by #2 = vr/2 — 9\ 
{fix = 35.3°, tan 2 6»i = 0.5, tan 2 6 2 = 2.0), with Am 2 = 7 x 10~ 5 eV 2 . The left hand panel 
shows the survival probability (averaged over 0.1 MeV energy bins) in vacuum for a source- 
detector separation of 180 km. The vacuum survival probabilities are identical for 6\ and 
#2- The right hand panel shows the daytime LMA oscillation solutions for solar neutrinos. 
Matter effects in the Sun cause the differences between the two solar solutions. 



gallium 0, [L5|, 0, Super-Kamiokande g7|, SNO CC pf, and SNO NC |J 



ex- 



periments, allowing the dominant p — p, 8 B, 7 B and CNO fluxes to vary freely but 
imposing the luminosity constraint (see ]5T] and references therein for details). We 
find that xLm solar = 53.9 for 2 dof (5 rates - 4 free fluxes + 1 constraint), which im- 
plies that non-standard solar model explanations of the solar data can be excluded 
at 7a. Adding to the analysis the evidence of antineutrino disappearance observed 
in the KamLAND experiment, we find that non-standard solar models invented to 
avoid solar neutrino oscillations are excluded by solar plus KamLAND experiments 
at 7.9a. 



3.1.2 Predictions for future experiments 

Table ||] presents the best-fit predictions and the expected ranges for some important 
solar neutrino observables. The results summarized in table |] correspond to the 
globally allowed oscillation regions shown in figure |j. We consider predictions for 
SNO, for a large SK-like neutrino-electron scattering experiment, and for a neutrino- 
electron scattering detector, BOREXINO JHJ or KamLAND p3|| , that observes the 
7 Be solar neutrinos. We also include predictions for a generic p — p neutrino detector. 
The notation used here is the same as in ref. [12j . 

The day-night asymmetry An-d is defined in terms of the Day and the Night 
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Observable 


b.f. ±1(7 


b.f ±3(7 


A N _ D (SNO CC) (%) 


q q+0.7 
°-°-0.6 


3.3111 


A N _ D (SK ES) (%) 


1.9 ±0.4 


i.9±?:I 


[R ( 7 Be)] v — e scattering 


0.64 ±0.02 


u - Dq: -0.04 


An n f 7 Be) f%) 


0+2-2 

u - u — 0.0 


0+^ 


[R ( 8 B)] v — e scattering 
(Tth = 3.5 MeV) 
(T^ = 5 MeV) 


0.46 ±0.03 
0.45 ± 0.03 


0.46 ± 0.04 
0.45 ± 0.04 


[p-p] v — e scattering 
(T th = 100 keV) 
(T th = 50 keV) 


0.697 ±0.023 
0.693 ±0.024 


0.6971852 
0.693±g52 



Table 2: Post-KamLAND Predictions. This table presents for future solar neutrino 
observables the best-fit predictions, and the la and 3a ranges, that were obtained from 
the global analysis shown in figure |||. 



event rates by the expression 

_ [Night - Day] 
An ~ D ~ 2 [Night + Day] " (3 ' 2) 

The prediction for the day-night asymmetry in the SNO CC measurement is one 
of the most interesting new results to come from the much smaller, post-KamLAND 
allowed regions (see discussion of earth matter effects in the LMA regime in refs. [[33 



S§). 

Within the unique la region (2 dof), the prediction is 

A N _ D (SNO CC) = 3.3+°; 7 %. (3.3) 

Within the 3cr region (2 dof), A N _ D (SNO CC) = 3.3tju% . 

An accurate day-night asymmetry measurement in the SNO CC mode could, in 
principle, help to discriminate between larger and smaller values of Am 2 . For values 
of Am 2 < 10~ 4 eV 2 within the currently allowed 3a (2 dof) region, 

v4 n _d(SNO CC) = 3.3+^%, (Am 2 < 10- 4 e\/ 2 ). (3.4) 
For values of Am 2 > 10 _4 eV 2 , the day- night asymmetry is small, 

A N _ D (SNO CC) = 0.6+5; 4 %, (Am 2 > 10~ 4 eV 2 ). (3.5) 
Unfortunately, it will be very difficult to reach the required experimental precision. 



13 



In table [| we also present the predicted day-night asymmetry at a SK-like 
neutrino-electron scattering detector with a threshold E t h = 5 MeV. For ES the 
day-night effect is decreased relative to pure CC interactions by the contribution of 
the neutral currents. 

The prediction for the reduced 7 Be v — e scattering rate, 

7 _ Observed v — e scattering rate 

^ ^ = BPOO predicted rate ' ( ' ' 

is remarkably precise and remarkably stable. Within the la allowed region, the 
pre-KamLAND prediction was [ 7 Be] = 0.64 ± 0.03, la |TJ. The post-KamLAND 
prediction for the currently allowed la region is 



'Be 



= 0.64 ±0.02. (3.7) 

post — KamL AND v ' 



The predicted event rate for the 7 Be rate experiment has a precision of ±3% in units 
of the rate calculated with the BPOO flux and assuming no neutrino oscillations. 

Now that LMA is the only oscillation solution allowed at 3a, the predicted 7 Be 
day-night asymmetry averaged over a year is too small to be observed. Within the 
current 3a allowed oscillation region, the limit on the predicted values of An-d is 

A N _ D ( 7 Be)| < 0.001. (3.8) 

We also include in table § predictions for the rate at which 8 B neutrinos will 
be scattered by electrons resulting in recoil electrons in BOREXINO or KamLAND 
above two different electron recoil kinetic-energy energy thresholds, Tth = 3.5 MeV 
and T t h = 5 MeV. The predictions are given in terms of the reduced event rate for 8 B 
neutrinos, [R( 8 B)], which is defined analogously to [ 7 Be] [see eq. ( |3.6|) 1 relative to the 
scattering rate predicted for the standard solar model flux with no oscillations. The 
predicted ranges of [R( 8 B)] are obtained from the allowed la or 3a regions in the 
three parameter space of Am 2 , tan 2 9, and /B,totai- Not surprisingly, the predicted 
8 B neutrino-electron scattering rates are essentially identical to the neutrino-electron 
scattering rates measured by the SNO |1| and Super-Kamiokande [[H]] experiments, 



which are 0.47±0.05 and 0.465±0.015, respectively. For plausible assumptions about 
how the BOREXINO detector will operate, one expects |~2| ~ 175 events above a 



3.5 MeV recoil electron energy threshold and ~ 110 events per year above a 5 MeV 
threshold. 

In addition, we have evaluated the predictions for a generic p—p neutrino-electron 
scattering detector. The reduced rate for this detector is defined by the relation 

r _ Observed v — e scattering rate 

^ " P * = BPOO predicted rate ' ^ ' ' 

We present the predicted rate for two plausible kinetic energy thresholds, 100 keV and 
50 keV. The predicted rate is precise and robust. The precision of the prediction is 
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±3%; the best-fit predictions have changed by less than 1% from the pre-KamLAND 
predictions [0 . 



The principal differences between the post-KamLAND predictions summarized 
in table [2] and the pre-KamLAND predictions given in table 2 of ref. [12j are a 
consequence of the disappearance of the LOW and the vacuum solutions. As a con- 
sequence, no significant day-night asymmetry nor seasonal variations are expected 
for 7 Be neutrino measurements by BOREXINO or KamLAND. For all other observ- 
ables, the main effect is the reduction of the la LMA ranges. Within the 3a LMA 
ranges, the changes are generally insignificant. 

3.2 Rates only analysis 

It is instructive to carry out a global analysis using only the total measured rates in 
the chlorine 0, gallium Super-Kamiokande |I7J, SNO CC [0, and SNO 



NC |]T^| experiments. In this simplified analysis, we ignore all spectral energy data 
and all measurements of the zenith-angle (or day-night) dependence. We let the total 
8 B neutrino flux be a free variable in the standard rates-only analysis. For simplicity, 
we consider pure active, or pure sterile, neutrino oscillations. Therefore, the rates- 
only analysis has three degrees of freedom (Am 2 , tan 2 8, and the 8 B neutrino flux) 
and five measured rates. 

The best-fit solution lies in the vacuum region. The oscillation parameters of 
the best fit solution are: Am 2 = 8 x Kr n eV 2 , tan 2 9 = 0.30(3.35) with xLn = 2 - 2 - 
The values of Xmm f° r the other well known solutions are: SMA (3.1) , LMA (3.8), 
and LOW (8.6), and pure sterile (23.1). All four solutions for active neutrinos, VAC, 
SMA, LMA, and LOW are allowed at 3a in the rates-only analysis. If only rates are 
considered, the VAC, SMA, and LMA solutions are essentially equivalent. 

We have done several tests to determine the robustness of this conclusion. The 
vacuum solutions always appear to be slightly preferred in the rates-only analyses, 
but the relative order of SMA-LMA can be reversed by slightly changes (such as 
using the original SNO CC data J3(| with a 6.75 MeV threshold rather than the 
more recent data |18|] with a 5.0 MeV threshold). 

This analysis illustrates the fact that the measurement of the essentially undis- 
torted recoil electron energy spectrum by Super-Kamiokande [171], which was con- 



firmed by SNO |T8|], was required in order to establish, prior to the KamLAND 



measurements, that LMA was the favored solution. Previous rates-only analyses 
have established the critical importance of the Super-Kamiokande spectrum analy- 
sis [0, |37j]. However, this is the first calculation that we know about in which all 
five of the currently available rate experiments have been included. 

3.3 KamLAND plus Solar Analysis 

It is convenient and conventional to discuss the experimentally determined 8 B neu- 
trino flux in units of the flux predicted by the standard solar model. Let /b, total 
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be the ratio of the total 8 B solar neutrino flux to the flux predicted by the BPOO 
standard solar model, i.e., 

/B ' t0tal = 0( 8 B) BPOO ' (3 ' 10) 



where 0( 8 B) BPO o = 5.05 x 10 6 cm'^- 1 |9 



We allow for the possibility that the 8 B neutrino flux that arrives at Earth 
contains a significant component of sterile neutrinos which we parameterize in terms 
of the active-sterile admixture parameter sin 2 77 (See Sec. |3.3.2|) . 

We calculate the global \ 2 by fitting to all the available data, solar plus reactor. 
Formally, the global x 2 can be written in the form 



4obai = xL ar (Am 2 ,tan^, /B>total ,sin 2 r ? ) + xL mL AN D (Am 2 , tan 2 



(3.11) 



There are four free parameters in Xgiobab Am 2 , tan 2 9, /b, total, and sin 2 r], although 
XxamLAND depends only on Am 2 , and tan 2 #. The most computationally intensive 
task in our analysis of solar plus KamLAND data, computing the multi-dimensional 
X 2 for all of the solar neutrino data, was calculated before the KamLAND announce- 
ment |J and stored. This data set amounted to about 9.4 gigabytes. We do not 



include the CHOOZ data [21] explicitly in eq. 3.11 since the combined solar data 
also excluded larger values of Am 2 . 

We have carried out a global analysis [12] of the combined solar neutrino data 



(80 measurements), see the discussion in section 3.3.1 below, and the KamLAND 
data (as described in section [2.2|) . We also consider the possibility that the 8 B neu- 
trino flux that arrives at Earth contains a significant component of sterile neutrinos. 
We describe in section |3.3.2| the aspects of the analysis that are related to sterile 
neutrinos. 



3.3.1 Solar analysis 

The details of our calculational methods are described in ref. [0. We treat the 
8 B solar neutrino flux as a free parameter, but adopt the predictions of the BPOO 
solar model |29| for the other neutrino fluxes and uncertainties. The free 8 B strategy 



adopted in the present paper was, in ref. |fL2| , the preferred analysis procedure among 
the three analysis strategies that were used earlier. 

We summarize briefly in this subsection the main features of our analysis of solar 
neutrino data. 

For the solar neutrino analysis, we use 80 data points. We include the two 



measured radiochemical rates, from the chlorine ][0| and the gallium [TJ|, 15, |I~6 | 
experiments, the 44 zenith-spectral energy bins of the electron neutrino scattering 
signal measured by the Super-Kamiokande collaboration ||17||, and the 34 day- night 
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spectral energy bins measured with the SNO [18|, [19], p8[ detector. The average rate 



for the gallium experiments [14, [TR [L6| has been updated to 70.8 ± 4.4 SNU by 



including the preliminary GNO data reported at Neutrino 2002 [|T| 

We take account of the BP00 |29|] predicted fluxes and uncertainties for all solar 
neutrino sources except for 8 B neutrinos. The Super-Kamiokande [|17j and SNO 



T9[] experiments can now determine the 8 B solar neutrino flux at a level of precision 
that is comparable with the precision of the BP00 prediction. Therefore, we treat 
the total 8 B solar neutrino flux as a free parameter to be determined by experiment 
and to be compared with solar model predictions. 

Following the SNO collaboration we have considered together (consistent 



with the way they are measured) the SNO charged current, neutral current, and 
electron scattering events. For each point in neutrino oscillation parameter space, we 
have computed the 34 day-night spectrum energy bins by convolving the computed 
survival probabilities with the appropriate neutrino cross sections, as well as the 
experimental energy resolution and sensitivity functions. 

We include the errors and their correlations for all of these observables. We 
construct an 80 x 80 covariance matrix that includes the effect of correlations be- 
tween the differ error as off diagonal elements. The errors and their correlations are 
described in in the Appendix of ref. and in section 5.4 of ref. [|1^ 



The analysis of the theoretical and experimental uncertainties is subtle. Many 
of the differences in the precise borders of the allowed oscillation regions that are 
reported by different groups may be traced to variations in the ways in which the 
errors are treated. We have performed the calculations accurately, although some 
of the refinements we have used are computationally demanding. As described in 
the Appendix of ref. [§], we include the energy dependencies and the correlations 
of the chlorine and gallium neutrino absorption cross sections. For both the Super- 
Kamiokande and the SNO experiments, we take account of the energy dependence 
and the correlation of the 8 B spectral energy shape errors, the absolute energy scale 
errors, and the energy resolution errors. 

3.3.2 Sterile component 

We have also considered the possibility that v e oscillates into a state that is a linear 
combination of active (v a ) and sterile (v s ) neutrino states, 

Ue^cosr/Ua + sin 77 z/ s , (3-12) 

where 77 is the parameter that describes the active-sterile admixture. This admixture 
arises in the framework of 4-z/ mixing . The total 8 B neutrino flux can be written 



<K 8 B) tot al = 0Oe) + 4>{Va) + </>(".), (3-13) 

where 0(V S ) = tan 2 rj x 0(z/ a ). 
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The solar neutrino data analysis depends on four parameters: the oscillation 
parameters Am 2 and tan 2 9, the total 8 B neutrino flux /b, total [ see ec i- ( ft.lOp l, and 
the sterile component characterized by sin 2 rj. We can also characterize the sterile 
contribution to the total flux in terms of /b, sterile, which is defined as 

0( 8 B, sterile) 

JB, sterile = T , gm ■ l^- 14 ; 

0( B)bP00 

We obtain from eqs. ( |3.13| ) and ( |3 . 1 2 ) the relation 

/B, s terUe = /b, total sin 2 r/ (l - P ee (E, Am 2 , tan 2 0, sin 2 77) J . (3.15) 

Thus for a given set of parameters Am 2 , tan 2 0, /b, total and sin 2 77, the inferred sterile 
flux parameter /b, sterile is an energy dependent, and time dependent, function. How- 
ever, the energy and time dependences are mild since, for the LMA solution, P ee is 
nearly constant for the 8 B neutrino energies that are accessible to direct measure- 
ment. 

In this paper, we account for the mild energy and time dependences by defining 
/b, sterile as the average sterile flux parameter inferred for the SNO experiment, i.e., 

iB.steriie = /b, total sin 2 r] (l - (P ee (E , Am 2 , tan 2 9, sin 2 r/)) S No) (3.16) 

where (P ee (E, Am 2 , tan 2 9, sin 2 t/))sno is 



tP (\ m ** m *ff\\ JdE » ^ SM ( 8B ' E u )a e {E u )P ee {E u , Am 2 , tan 2 Oj 
\^ee\^rn ,tan t/j; SNO = Bssm . {o.ifj 

Tic 



L SNO 

Here P ee is the average survival probability and -Rg^o is the predicted SSM |29| rate 
for the CC SNO measurement |I§L pOl in the absence of neutrino oscillations, 



Rl™ = J dE u SSM ( 8 B, E v )a e {E v ). (3.18) 



The cross section a e (E v ) ^0], ^T| is the weighted average cross-section for the CC 
reaction, v e + 2 H — > e + + n + p, for the neutrino energy E v . The average includes 
the experimental energy resolution function, Res(T, T'), where T(T') is the measured 
(true) recoil kinetic energy of the electron. Thus 

a e (E u ) = dT dT' Res(T, T') j . (3.19) 

JT min Jo dT' 

The lower limit, T m i n , in the integral in eq. Q3.19 ) is taken here to be the threshold, 
5 MeV, used by the SNO Collaboration in ref . |18| . The calculated value for the CC 
rate is not sensitive to the assumed value of T max , as long as T max > 17 MeV. 

We construct Xsoiar(^ m2 > tan 2 9, /b, total, sin 2 r/) in four-dimensional neutrino pa- 
rameter space. The results shown in figs. |3] and f| are the two dimensional allowed 
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regions in Am 2 , tan 2 9 that are obtained after marginalizing xLiar an d Xgiobai [de- 
fined in eq. ( |3.11|) 1 over /b, total an d sin 2 ?]. The allowed domain of the oscillation 
parameters at a given CL is defined from the condition 

xLrginaiizcdlA^^an 2 ^) < X Ln + A% 2 (2, CL) (3.20) 

where Xmm = 81-2 is the global minimum for the solar plus KamLAND analysis 
(for 89 = 80 + 13 — 4 dof) in the four-dimensional neutrino parameter space and 
Xmarginaiized(^ m2 5 tan 2 9) is obtained by minimizing with respect to /b, total and sin 2 // 
either x 2 olar or Xgiobai f° r eacn value of Am 2 and tan 2 9. 

Conversely, in section [| (and |) we determine /b, total (and the sterile contribu- 
tion) by minimizing Xgiobai with respect to Am 2 , tan 2 9, and sin 2 ?] (or /b, total)- The 
allowed ranges for each of these parameters at a given CL are obtained from the cor- 
responding condition for 1 dof. For instance, the allowed range of /b, total is obtained 
by applying condition eq. ( |4.1|) , and similarly for sin 2 ?? (or /b, sterile)- 



4. The total 8 B solar neutrino flux 

One can determine the allowed range of the total 8 B neutrino flux using the results 
of the KamLAND reactor neutrino experiment and the results of solar neutrino 
experiments ref. . The conceptually simplest approach is to use just the KamLAND 
determination of the neutrino oscillation parameters (cf. figure [I]) together with the 



measured CC event rate in the SNO |18| experiment to infer the total B solar 



neutrino flux. However, this KamLAND-only approach, while very intuitive, does not 
make use of the greatly reduced allowed regions obtained when solar and KamLAND 
data are combined (cf. figure and figure f|). With the present accuracy of the 
KamLAND data set, the KamLAND-only determination yields an imprecise value 
for the total 8 B neutrino flux (la uncertainty ~ 30%). 

Therefore, we present in this section the more accurate total 8 B solar neutrino 
flux that is derived from the four-dimensional global analysis of solar plus KamLAND 
neutrino data (see definition of Xgiobai i n ec L- (|3-11|) and discussion in section |3~3| ). Our 



results for /b, total are summarized in columns three and four of table |3|. 

In determining /b, total, we can use all of the available solar neutrino data as 
well as the KamLAND measurements. This method is much more computationally 
demanding than the more intuitive KamLAND-only approach since in the global 
approach one includes, in addition to the KamLAND data, all of the 80 solar mea- 
surements and their associated (sometimes correlated) uncertainties. Formally, we 
calculate (as in section ||) a global x 2 , see the definition of Xgiobai defined in eq. (|3.11| ). 

We obtain the allowed range of /b, total, the total 8 B solar neutrino flux in units 
of the BP00 predicted flux [see definition in eq. ( |3.10|) 1 by marginalizing Xgiobai with 

global 



respect to Am 2 , tan 2 9, and sin 2 ?]. To be explicit, we minimize Xgiobai f° r eacn value 
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Am 2 


tan 2 6> 


x 2 


/b, total 


3<7 /b .total 


,/b, sterile 


3<7 /B,sterile 


7.1 x 10" 5 


4.5 x KT 1 


81.2 


1.00 


[0.81,1.36] 


0.0 


[0.0,0.45] 


1.5 x 1(T 4 


4.3 x 10" 1 


87.3 


0.88 


[0.83,0.92] 


0.0 


[0.0,0.08] 



Table 3: The total and the sterile 8 B solar neutrino flux. For the two currently 
allowed LMA islands in the space of neutrino oscillation parameters that are shown in 
figure [|, the table presents the best-fit values for the total 8 B solar neutrino flux /b (total) 
and the sterile component of the 8 B neutrino flux, /b (sterile). The values are given in 
units of the predicted standard solar model 8 B neutrino flux ||29[|. We also present the 
associated 3a allowed ranges for the global solar plus KamLAND analysis. The different 
allowed regions are labeled by the best-fit point in oscillation parameter space. 



of /b, total with respect to Am 2 , tan 2 9, and sin 2 rj. The corresponding range of /b, total 
at a given CL is obtained from the condition 

^marginalized (/b, total) < xLn + Ax 2 (l,CL). (4.1) 

We find, after extensive numerical explorations of the four-dimensional space of 
Xgiobai' that the current best empirical estimate for /b, total is 

/ B ,total = 1.00(1 ±0.06) (la). (4.2) 



The result given in eq. ( |4.2| ) agrees with the standard solar model prediction p9 



In table |3], we show the ranges for /b, total for within each of the LMA islands in 
figure 0L 



5. The sterile component of the 8 B solar neutrino flux 

The sterile component of the 8 B solar neutrino flux can be determined directly by 
subtracting the active 8 B neutrino flux, which is given by the SNO NC measure- 
ment [rS , from the total 8 B neutrino flux determined in section |] (cf. refs. 0, [|). 
This subtraction can be made in a model independent fashion using the oscillation 
parameters obtained from the analysis of the KamLAND data alone. However, this 
direct method, while conceptually simple, yields an imprecise determination of the 
sterile component of the 8 B neutrino flux (lex uncertainty of order 18%). 

The more powerful approach is to include all of the solar plus KamLAND 
measurements (cf. section |j). In this approach, we first evaluate the global \ 2 i 



Xgi OD ai(Am , tan /b, total, sin 77), which is denned in eq. (|3.11|) . Not only are all 
of the relevant data used in calculating Xgiobab Du t a lso this method properly in- 
cludes the correlated uncertainties. Thus calculating Xgiobai is the preferred statistical 
method as well as the more powerful (but computationally intensive) method. 
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Figure 6: Limit on active-sterile admixture. The figure shows the active-sterile 
admixture sin 2 rj as a function of /b, total, the total 8 B flux relative to the BPOO predicted 
standard 8 B neutrino flux . The contours shown in figure || are 90%, 95%, 99%, and 99.73% 
(3cr ) confidence limits. 



As discussed in refs. @, 20|, the allowed sterile component in the 8 B neutrino 
flux is strongly correlated with the allowed departure from the standard solar model 
prediction of the total 8 B neutrino flux. We illustrate this correlation in figure || 

We show in figure |6| the allowed regions in the /b, total — sin 2 ?] parameter space 
after marginalizing with respect to Am 2 and tan 2 9. For the sake of clarity, we show 
separately the results of marginalizing for both of the LM A islands shown in figure f|. 

In figure [7], we display the corresponding allowed regions in terms of /b, total — 
/b, sterile (after marginalizing with respect to Am 2 and tan 2 #). 

Finally, marginalizing with respect to Am 2 , tan 2 8, and/s, total, we obtain the 
allowed range of the active-sterile admixture 



sin 2 r] < 0.13 (0.52) 



(5.1) 



at la (3a). 

Correspondingly for /b, sterile, we show in columns five and six of table |3| the best- 
fit values and the 3a ranges for the local minima of the two currently allowed LMA 
islands. The la range from this global analysis is 



B. sterile 



0.00+™ 



(5.2) 



We use eq. ( |5.2| ) as our best-estimate of the sterile contribution to the observed B 
solar neutrino flux. 
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Figure 7: The Sterile- Total flux correlation . The figure shows the sterile component, 
/b, sterile as a function of /b, total; the total 8 B flux, both relative to the BPOO predicted 
standard 8 B neutrino flux . The contours shown in figure [7| are 90%, 95%, 99%, and 99.73% 
(3cr ) confidence limits. 

Is the limit on / B , sterile given in eq. ( |5.2| ) (or on sin 2 rj in eq. ( |5~T| ) changed signif- 
icantly if a more general mixing of sterile neutrinos is considered? We have investi- 
gated numerically cases in which a second Am 2 affects solar neutrino oscillations. We 
have adopted as a paradigm the model discussed in ref . |§ , in which there is a third 
neutrino mass scale, Am 2 ~ 10 _6 eV 2 , related to a sterile neutrino. We find that, 
within the accuracy indicated by eq. ( |5.2|) , our best-estimate of the sterile contri- 
bution (and the associated uncertainty) of the 8 B solar neutrino flux is not affected 
by the presence of this lower mass scale. The main effect of this low-mass sterile 
neutrino scenario is that the sterile neutrinos add an extra energy dependence to the 
survival probability and to the sterile fraction. So the limit given in eq. ( |5.2| ) for the 
8 B solar neutrino flux will not apply to the sterile component of the lower energy 
p — p and 7 Be solar neutrino fluxes. The sterile component for the p — p and 7 Be 
neutrino fluxes must be determined by a combination of CC and v — e scattering 
experiments at low energies (< 1 MeV). 

6. Summary 

The allowed oscillation regions that are determined by including the initial Kam- 
LAND results together with the previously known solar neutrino results (see 
figure ^) populate 6 orders of magnitude less area in neutrino oscillation parameter 
space (Am 2 , tan 2 9) than the allowed oscillation regions that are determined with 
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just solar neutrino data (see figure 0). The LMA solution is now the only allowed so- 
lution at 3a. The LOW, Vacuum, and SMA solutions are disfavored at, respectively, 
4.7<T, 4.8cr, and 6.1<r. All alternative explanations to the solar neutrino problem that 
do not also predict a deficit in KamLAND are now excluded at 3.6cr. Combining 
solar neutrino data with the KamLAND measurements, non-standard solar model 
explanations of the solar neutrino problem that do not include new particle physics 
(e.g., neutrino oscillations) are excluded at the 7.9a CL. Maximal mixing, tan 2 9 = 1, 
is disfavored at 3.1a. The results in this paragraph are described in sec [3.1.1| . 

At the present level of accuracy, matter effects are not important in determining 
the allowed oscillation regions permitted by the KamLAND experiment. However, 
when the precision of the experiment is improved, matter effects must be included 
(see sec. |2.3| ) . Within the currently allowed oscillation regions, the maximum change 
is 4% in the calculated event rate for the KamLAND experiment for a specific value 
of Am 2 and tan 2 9. 

The LMA prediction of the v — e scattering rate for 7 Be solar neutrinos is precise, 
0.64±0.03 within the la (2 dof)currently allowed oscillation region (see section |3.1.2| ), 
in units of the rate predicted by the standard solar model and assuming no flavor 
conversion. The estimated uncertainty in the predicted standard solar model 7 Be 



solar neutrino flux, ±10% p9 |, is much larger than the current uncertainty in the 
predicted scattering rate due to neutrino oscillation parameters. Thus BOREXINO 
and KamLAND can provide a crucial test of solar model predictions by measuring 
accurately the 7 Be v — e scattering rate. 

Table of section [3 . 1 . 2| summarizes the predicted rates for v— e scattering by p— p 
as well as 7 Be ( and 8 B) neutrinos. Within the la (2 dof) allowed region, the predicted 
p—p scattering rate is ~ ±2% due to imperfectly known oscillation parameters, which 
is about twice the currently estimated uncertainty in the solar model prediction of 
the p — p neutrino flux. If our current ideas regarding the accuracy of the solar model 
predictions are valid and if the popular neutrino oscillation framework is correct, then 
a measurement of the p — p solar neutrino will provide a precise determination of the 
mixing angle 9. 

It will be difficult with the existing SNO detector jn| [Tj| to measure the pre- 
dicted SNO CC day-night asymmetry: A N _ D (SNO CC) = 3.3t^%( la) (see 
eq. |3.3p . 

The event rate measured at KamLAND, (eq. |2.1| ) lies within the 1-sigma pre- 
KamLAND expectations based on purely solar neutrino experiments (eq. |2.2| ). This 
excellent agreement indicates that the main features of solar neutrinos are beginning 
to be well understood. 

The most accurate method to determine the total 8 B flux is to combine together 
all of the solar as well as KamLAND experimental data. Using this computationally 
intensive method, we find (see section ^) that the total 8 B solar neutrino flux is 
1.00 (1 ± 0.06) in units of the standard solar model predicted flux |29| . 
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The sterile component of the 8 B neutrino flux can be determined similarly by 
using all of the available solar and KamLAND data. The allowed range of the active- 
sterile admixture is sin 2 77 < 0.13. Correspondingly, our best estimate for the sterile 
component of the 8 B solar neutrino flux is O.OOloioo m units of the standard solar 
model predicted flux (see section H). 

Our results for the total 8 B solar neutrino flux, and for the sterile component of 
this flux, are summarized in table ^. The present determination of the total 8 B solar 
neutrino flux reduces the la experimental uncertainty in this quantity by a factor of 
two with respect to the previous most accurate determination [0 . The experimental 
bounds on the total 8 B neutrino flux are now more than a factor of two smaller that 



the uncertainties in the predicted standard solar model flux ||29|| . This is the first 
ocassion on which the experimental errors on a solar neutrino flux are smaller than 
the solar model uncertainties. The uncertainty in the sterile neutrino component of 
the 8 B neutrino flux represents a factor of two improvement over the previous most 
stringent bound. 

We conclude with an important methodological question. Do the results given in 
this paper depend in a significant way upon details of the analysis of the KamLAND 
data? This is a natural question to ask since the KamLAND results represent a 
crucial but low statistics addition to the solar neutrino data. Fortunately, different 
analysis strategies lead to essentially the same quantitative inferences regarding solar 
neutrino properties and predictions. 

We have studied the robustness of our quantitative conclusions to different plau- 
sible variations in the statistical analysis of the KamLAND data. Most of these 
variations made no discernible difference in the final answers that were obtained 
from the global solar plus KamLAND analysis. For quantities such as the total 8 B 
neutrino flux, the sterile component of the 8 B neutrino flux, and the predictions for 
observables in future experiment, there were no significant differences. For speci- 
ficity, we report here the results of a x 2 analysis of the KamLAND data which led 
to the largest changes of any of the alternative analysis strategies we adopted. In 
this x 2 analysis, we used the first 5 energy bins above 2.6 MeV visible energy as pre- 
sented by the KamLAND collaboration but collected all the events above the fifth 
bin (E > 4.725 MeV) into a single bin (or into two bins) in order to have enough 
events in each bin to assume Gaussian errors. We found that the allowed regions of 
the KamLAND-only analysis are slightly modified, but in ways that would be noticed 
only by the trained and expert eye (not to the casual reader). However, the effect of 
using a standard x 2 with Gaussian errors is even less important for determining the 
global, solar plus KamLAND, allowed regions. In particular, this alternative analysis 
changes /b, total by ±0.01 (la) and /b, sterile by ±0.005 (la). The maximum change 
in the range of predicted solar neutrino observables in table [2] occurs for the _R( 8 B) 
and amounts to 2% at la in the largest allowed value. For all other observables, the 
changes are less than 2% in the 3a range. 
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